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I. INTRODUCTION 


A. BACKGROUND 


The complex nature of many common heat transfer situations causes most of the 
work in thermal fluids to deal with some portion of experimentation. In more complex 
(and typically, more realistic) systems, analytical solutions are not possible. Whether it is 
conducting the actual experiments, reviewing the results, using the results, or correlating 
the results, it is not uncommon for undergraduate students all the way up to corporate 
engineers to be dealing with experimentation of thermal fluid systems. 

With more and more complexity in fluid systems and experimentation, correlation 
is becoming more and more difficult. Most experiments involve varying fluid 
properties, variable heat transfer conditions and a host of geometric configurations. With 
all of these parameters, concise methods of graphical display are impossible and 
correlations are becoming far too complex for use without a computer program. Because 
of the assumptions made by these correlations (such as constant fluid properties), 
inaccuracies and error occur. Even the simplest, most limiting case, correlation 
development is quite an undertaking and often leads to equations that must be evaluated 
using computer programs. 

In complex situations, the use of traditional, exponential forms to correlate the 
experimental data results in predictive expressions with inherently large uncertainties. 
Using the neural network technique, this raw data becomes easier to use as a predictive 


and design tool. 


Neural networks provide an alternative to both graphically represented data and 
the use of complex correlation equations. The neural network technique provides distinct 
advantages over other approaches as well as the highest degree of accuracy of the three 
methods discussed (graphical displays, mathematical correlations, and neural networks). 
Although the neural network displays usefulness in control, design, and prediction, only 


the area of prediction of heat transfer parameters will be discussed. 


B. THE ARTIFICIAL NEURAL NETWORK 


Though it is not the focus of this study to provide a detailed discussion of neural 
networks, the network that is used will be described so that its use will be understood and 
enough references will be given if more detail is desired. 

The neural network (and loosely artificial intelligence) has been defined in many 
ways. For this study, the definition provided by Sage is most appropriate, "the 
development of paradigms or algorithms that require machines to perform cognitive 
tasks, at which humans are currently better" [Ref. 4]. It may seem inconsistent to use this 
definition while stating, as above, that the neural network in this study is being used to 
replace correlations: replacing the computer (rather than the human) with the neural 
network. However, the neural network will tackle groups of similar correlations. In 
essence, the neural network will be replacing the human selection of the correlation and 
the correlation itself. 

While it has flourished in other disciplines, the thermal sciences have been a bit 
reluctant to accept the neural network technique as a method of prediction [Ref. 3}. This 


is most likely due to the lack of system analysis necessary in order to develop a solution. 





Like the human brain, the neural network has no need for an equation in order to 
develop a solution. This can be best understood through an example from M. Sen of the 
University of Notre Dame. Take, for example, your car. When you turn the wheel (an 
input), you do not know how many degrees your front tires will turn. You do, however, 
know how your car will react when rounding the turn (the output). From human 
standpoint, it is unnecessary to know the degree value that your tires turned, but only the 
end result of rounding the corner. For a heat exchanger, this would translate as only 
having the hot and cold fluid temperatures at the inlet and outputting the output 
temperatures--no other information is necessary. This is much different than any 
analytical approach and thus, the community is reluctant to invoke the neural network 


technique. 


C. PREVIOUS WORK 


This study, though not a direct continuation of their work, has been inspired by 
the work done by Sen and Yang of the University of Notre Dame. They have had success 
in the simulation and control of heat exchanger performance using neural networks. 
Other heat transfer applications of the neural network technique have been conducted. 
More specifically, in 1996, Jamunathan et al. [Ref. 9] worked in the area of liquid crystal 
thermography. In 1994 and 1995, Lavric et al. [Ref. 10] studied the design of a finned 
tube heat exchanger. In 1994, Huang and Nelson [Ref. 11] determined the delay time for 
HVAC plants and in 1990, Ding and Wong [Ref. 12] attempt the control of a simulated 
hydronic system [Ref. 3]. The use of neural networks in thermal fluid problems is a rather 
"young" subject as shown by the aforementioned exhaustive list of known studies. 


However, the Artificial Neural Network technique (ANN) is a well-established technique 


[Ref. 4]. This study breaks away from other neural network/thermal fluid studies by 


displaying its usefulness stnctly in any area involving the correlation of data. 


D. OBJECTIVE 


The objective of this study is to investigate the feasibility of employing the neural 
network technique as a method of using experimental data to predict heat transfer 
behavior. Currently, data 1s acquired by experimentation, collected en mass, and then 
correlated to one or more of the controllable inputs using some fluid mechanics insight 
and some mathematical know-how. Experimental uncertainties in the data accumulation 
are coupled with the inherent uncertainties in the mathematical correlation. The 
correlation, in order to be useful, 1s typically called upon to accommodate a range of 
inputs that leads to more error. Incropera and DeWitt [Ref. 2] comment that the 
correlations should not be viewed as sacrosanct, "Each correlation is reasonable over a 
certain range of conditions, but for most engineering calculations one should not expect 


accuracy too much better than 20%." 


It is the goal, then, to make the predictions of thermal fluid behavior more reliable 
(i.e. less error), less reliant on assumptions, and provide easier methods of evaluating 
these predictions. With neural networks all of the above goals are realized. However, as 
an additional goal, this report attempts to familiarize the thermal fluid community with 


the unique capabilities of ANN as a predictive technique. 





II. THE NEURAL NETWORK 


A. GENERAL NEURAL NETWORKS 


As mentioned previously, this report does not intend to be a survey in neural 
networks. However, it will offer a brief explanation and will make reference to many 
other sources of such information. For a mathematical and historical background of 
neural networks reference the text by Haykin [Ref. 4]. The paper by Sen and Yang [Ref. 
3] provides an extensive review of the applications of artificial neural networks and 
genetic algorithms in thermal engineering. 

Although many different types of neural networks exist, this study and most 
engineering applications primarily invoke the use of the supervised fully connected single 
and multi-layer feedforward network with backpropagation. As its name implies, neural 
networks are based on the biological make-up that runs the brain, whether it is human or 
non-human. Like the brain, the neural network learns, recognizes information, and 
makes generalizations on the input it receives [Ref. 3]. The essence of the method 1s to 
"train the network" with inputs that correspond to known outputs. Through a process that 
is roughly labeled "experience," the neural network 1s able develop the understanding of 


what output results from certain inputs. 





Output 
Input First Output 
Layer Hidden Layer 
Layer 
Figure 1: Schematic of a General Backpropogation Network. 
Ie General Layout 


Figure 1 presents the layout of a single hidden layer fully connected feedforward 
network. This layout represents only one type of network, albeit the most important 
network for this study. Each gray circle represents a neuron of the network and 
implements a function to determine the output of that neuron. The input layer shows that 
this network has four inputs, one hidden layer (this layer is neither an input or output 
layer), and an Output layer with a set of neurons. The hidden layer and output layer 


invoke a function that determines each layer's output. 


6 


Input Neuron without bias 


a= f(wp) 


Figure 2: Schematic of Neuron without Bias. 


pm Neuron 


A neuron with a single scalar input and no bias is shown in Fig. 2 and with bias in 
Fig. 3. The scalar input p is multiplied by the weight function w and remains a scalar 
value wp. For the neuron without bias, the wp is the only argument for the transfer 
function. For the neuron with bias, the argument for the transfer function 1s altered by 0d. 
The transfer function £ typically a sigmoid or a linear transfer function (discussed in the 
following section), is set before the network is trained. It takes the argument wp (or 
wp+b) and produces the output, a. Note that w and 6 are both adjustable. These are the 
quantities that are adjusted so that the network can produce output that is consistent with 


the training [Ref. 7]. 


Input Neuron with bias 





a= f(wp +5) 


Figure 3: Schematic of Neuron with Bias. 


Se Transfer Functions 


Each neuron has an input/output characteristic and implements a transfer function. 
The two transfer functions used in this study are the linear transfer function (shown in 
Fig. 4) and the log-sigmoid transfer function (shown in Fig. 5). The linear transfer 
function calculates the neuron's output by returning the value passed to it. The log- 
sigmoid transfer function takes the input, which is any real value and returns with an 
output between the range of 0 and 1. The algebraic form is given in Eq. 2.1. 


] 
l+e” 





f(n)= (2.1) 


Networks with biases, a sigmoid layer, and a linear output layer are capable of 


approximating any function with a finite number of discontinuities [Ref. 7]. 
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Figure 4: Linear Transfer Function. 





Figure 5: Log-Sigmoid Transfer Function. 


Input Layer of Neurons 





a= f(Wp+b) 


Figure 6: A Detailed Schematic of One Layer of a Network. 


Figure 6 shows a single layer network of S log-sigmoid neurons having A inputs. 
In order to achieve an output outside of the range of 0 and 1, another of layer of linear 
transfer functions would be necessary. For this type of network the input of the next 
layer would be the a vector. This type of network is shown in Fig. 7. A vector of input 


p is the argument of the hidden layer transfer function. The output of this layer, a, is the 
argument for the output layer transfer function. The vector a,is the output of the second 


layer and the network altogether. 
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Layout of Input, Hidden and Output Layers [Re 


Figure 7 


1] 


Also in Fig. 7, notice the transfer functions are given below the schematic of the 
network. The hidden layer transfer function is a log-sigmoid and the output layer is the 


linear transfer function. Each layer has a bias, b, and a layer (or input) weight LW (or 


IW) [Ref. 7]. 


4, Training 


The training of the network requires a set of examples of proper network behavior 
from which the network can make proper pattern associations. The set of examples has 
both the vector input (or matrix of inputs) and the output corresponding to the input (the 
targets). During the training, weights and biases of the network are iteratively adjusted to 
minimize a function set as the measure of network performance. Typically, the gradient 
of the performance function is used to determine the adjustments of the weights and 
biases. The gradient is determined using a technique which involves performing 
computations backwards through the network. This backpropagation invokes 
optimization techniques (such as Newton methods, secant method or conjugate gradient) 
to compute the gradients for multi-layered networks. Input vectors and the corresponding 
Output vectors are used to train a network until the approximation of a function is 
reached. Properly trained backpropagation networks typically approximate functions 
well even when presented with inputs different from the training data [Ref. 7]. This 
generalization allows the training of a network to be done with a finite set of data, rather 


than an exhaustive data bank. 


B. THE SPECIFIC NEURAL NETWORK 


In this section, the primary network that was used in the study will be discussed. 
Like the general network described above, this network has the same basic architecture: 
an input layer, a hidden layer with a set number of neurons, and an output layer. It 
includes both a log-sigmoid (hidden layer) and a linear transfer function (output layer), 
and was trained using a set of example data. This type of network is shown in Fig. 7. In 


order to make the network more efficient the training data was preprocessed. 


Ae Preprocessing 


A neural network's training can be made more efficient if certain preprocessing 
steps are performed on the network inputs and targets. The training data inputs and 
targets were preprocessed by normalizing them by their respective mean and standard 


deviation. 


ya Initialization 


Once the network architecture is set, the network variables (weights and biases) 
are initialized. This gives values to the weights and biases so that the first 
backpropagation iteration, or epoch, can begin. Typically one epoch of training is defined 
as a single presentation of all input vectors to the network. The network is then updated 


according to the results of all those presentations. 


3: Simulation 


With the network initialized, a simulation can be done. This will confirm that the 


network is setup to receive the input, apply the weights and biases and produce an output. 


However, without iterating to find the optimal weights and biases, the output of the 


network will not approximate the targets with any accuracy. 


4. Training 


With the assumption that the network was setup with enough neurons in the 
hidden layer, the training will be able to adjust the weights and bias in order to 
approximate any function with a finite number of discontinuities. The backpropagation 
algorithm updates the weights and biases in the direction in which the performance 


function decreases most rapidly. The iteration can be written [{Ref. 7] 
Xp.) =X, ~ A,B, (2) 
where x, 1S a vector of current weights and biases, g, is the current gradient, and a, is 


the learning rate. The training is then done in batch mode where all of the inputs are 
applied to the network before the weights are updated. The gradients calculated at each 
training example are summed to determine the adjustment of the weights and biases. 

In order to find the amount to adjust the weights and biases the Levenberg- 
Marquardt algorithm is used. By employing this technique, a second-order training speed 
is achieved by approximating the Hessian matrix. The Hessian matrix is vital to 
backpropagation learning dynamics. According to Haykin [Ref. 4]: 


The Hessian Matrix of the cost function €,(w), denoted by H its 
defined as the second derivative of €_(w) with respect to the 
weight vector w, as shown by 


Ate o"So(W) 


~ 


The Hessian matrix plays an important role in the study of neural 
networks; specifically, we may mention the following: 


l. The eigenvalues of the Hessian matrix have a profound 
influence on the dynamics of back-propagation learning. 

pi The inverse of the Hessian matrix provides a basis for 
pruning (i.e. deleting) insignificant synaptic weights from a 
multi-layer perceptron, as discussed in Section 4.15. [in 
Haykin] 

3. The Hessian matrix is basic to the formulation of second- 
order optimization methods as an alternative to back- 
propagation learning, as discussed in Section 4.18. [in 
Haykin] 


The Levenberg-Marquardt algorithm approximates the Hessian matrix in the following 


iteration [Ref. 7]: 


x,., =x,-[J7J +1] J7e (2.2) 


and 


where e is the vector of network errors with respect to the weights and biases, H is the 
Hessian matrix, J is the Jacobian matrix of the weights and biases, 1 representative value 


of the performance function. 
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Il. CORRELATING THE GRIMISON DATA 


A. DESCRIPTION OF EXPERIMENT 


To show the usefulness of neural networks in thermal fluids applications, a set of 
training datais needed. Although there are no real requirements on what type of data is 
desired, a few guidelines helped the selection of the data set. For this experiment it was 
desired that the data set be well-known, easily (and relatively quickly) entered ona 
computer, and complicated enough to display the neural network's usefulness. For these 
reasons the Grimison data [Ref. 8] was chosen. 

Although the Grimison data are not numerically accurate for each arrangement, 
the network was trained as if the data were accurate (i.e. the data was used as the training 
data). Gnmison's relations were developed from the original experimental investigations 
reported by Huge [Ref. 6] and Pierson [Ref. 5]. The data by Huge and Pierson would be 
ideal for use as the training data. However, in each of their reports, data 1s recorded 
graphically which makes accurate use of their data quite difficult. For this reason the 
Grimison correlation data is used instead. 

Grimison has obtained a correlation of the form [Ref. 2] 

N, 210 


Nun = © Renee 2000 < Rep na, < 40,000 (3.1) 
heen 7 


where C, and are listed in Table 1. 





Figure 8: Schematic of Aligned Configuration. 


The tube arrangements are shown in Figs. 8 (aligned) and 9 (staggered). The 


Reynolds number is given by 
Nene = (3.2) 


where p is the fluid density and w is the dynamic viscosity. V_.. for each tube 


arrangement is given by [Ref. 2] 


a _ ay, (for aligned, see Fig. 8) (3.3) 


max 
or for the staggered configuration, the maximum velocity may occur at the transverse 


plane A, or the diagonal plane A, (see Fig. 9). It will occur at A, if 


NS. IO (Se 18), (3.4) 


18 


and the value will be given by [Ref. 2] 


If V_,. occurs at A,, it may be computed like the aligned configuration in Eq. 3.3. 


Voi oat 7, ee 
2(S>p =i) 


Sr 


V (for staggered, see Fig. 9) 


(3.5) 


Table 1: Table of Values to be Used with the Grimison Correlation [Ref. 8]. 


See 
Aligned 
eS 
1.50 
2.00 
3.00 
Staggered 
Nhe) 
1.50 
2.00 
3.00 


C, 


0.348 
O96) 
0.418 
0.290 


0.518 
0.451 
0.404 
0.510 


eZ 


m 


W592 
0.586 
Fo 0 
0.601 


0.556 
0.568 
roe 
0.592 


Obras 
0.250 
209 
Weg] 


0.505 
0.460 
0.416 
0.356 


0.608 
0.620 
0.602 
0.584 


0.554 
0.562 
0.568 
0.580 


0.100 
0.101 
e229 
0.374 


O19 
0.452 
0.482 
0.440 


0.704 
0.702 
0.632 
0.58] 


0.556 
0.568 
0.556 
0.562 


0.0635 
0.0678 
0.198 
0.286 


O22 
0.488 
0.449 
0.428 


O52 
0.744 
0.648 
0.608 


0.562 
0.568 
O70 
0.574 


Table 1 and Grimison's associated equation (Eq. 3.1) show that this data set has 


four inputs: the transverse spacing S,, the longitudinal spacing S, , the Reynolds number 


Re and whether the tube rows of a bank are either staggered or aligned and one output: 


the Nusselt number Nu, . Four values different values were used for S, and S, (See 


Table 1). Re was varied between 2,000 and 40,000 at increments of 2,000 for use in Eq. 


3.1. Each of these data sets was used for both the staggered and the aligned configuration 


for a total of 640 training data sets. 





Figure 9: Schematic of Staggered Configuration. 


B. PRESENTATION AND USE OF DATA 


Using the arrangements shown schematically in Figs. 8 and 9, Pierson and Huge 
independently developed sets of data of Nusselt number as a function of Reynolds 
number for different arrangements. As described their reports, much effort was made to 
maintain accuracy throughout the experimentation. Their data, because quite valuable at 
the time, was reported exclusively in graphical form (See Fig. 10). Trends are 
immediately apparent, however. Nusselt numbers were shown to increase "roughly as 


the two-thirds power of Reynolds number" [Ref. 5]. Yet no correlation was made. 
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Figure 10: Example of How Pierson Presented Data in Graphical Form [Ref. 5]. 


Although the data was extremely useful in its own night, not until Grimison made a 
correlation did the data find usefulness in the investigation of tube banks. With the use of 
the correlation (and the associated table of C, and m values, See Table 1), the Nusselt 
number could be predicted without experiment or interpolation of Pierson's or Huge's 
data. Grimison, however, admittedly modified the original observed and measured data. 
Even after these alterations, the bulk of the data was reported to be on the order of +5% 

with the staggered banks "somewhat greater" [Ref. 8]. 

Using the Grimison data as the target (for reasons described above), the neural 
network is able to provide much more accurate results. The network is trained with one 
set of data (Nu = f(Re,S,,5,,config)). In order to test its ability to predict system’s 


behavior, a different set of input is given to the network and the output of the network 1s 


compared to the output of the correlation. 
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Ile Reynolds Number 
Figures 11 and 12 show that the network could make accurate predictions given 
Reynolds numbers different from those with which it was trained. Since there are four 


inputs, the plot displays the output of the neural network versus what the Grimison data 
would provide. A perfect correlation would produce a 1:1 ratio (a 45° line). In the worst 


case, 


%Error = te =O (me) (3.6) 


In Tables 2 and 3, the aligned and staggered % error is shown. Each value corresponds 
to the average over the range of Reynolds number (2000 < Re < 40,000 ). 

The correlation was reached relatively easily and quickly after developing a 
suitable network to train. The hidden layer consisted of 200 neurons and one neuron in 


the output layer. 
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Figure 11: Plot of Neural Network Output (ANN) vs. the Values Attained by the 
Grimison Correlation. 
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Figure 12: Plot of Neural Network Output (ANN) vs. the Values Attained by the 
Grimison Correlation. This is a Portion of the Previous Figure Shown in More 
Detail. 


Table 2: Average % Error (Over All Re) for Every Configuration of the Aligned 
Configuration. 


2S 1.50 2.00 3.00 





0.106 0.0791 0.0969 0.0717 
0.107 0.0418 0.0738 


2.00 | 0.0826 0.0721 0.0727 0.0776 
3.00 0.106 0.604 0.0630 0.179 











in 





Table 3: Average % Error (Over All Re) for Every Configuration of the Staggered 
Configuration. 


Pet 


Ss ae or Vr 
Ps 0.0426 {| 0.0512 | 0.0228 | 0.0196 










2 Longitudinal Spacing (S, ) 
Figure 13 shows that the network could make accurate predictions given 
longitudinal spacing (S,) different from that with which it was trained. Since there are 


four inputs, the plot displays the output of the neural network versus what the Grimison 
data would provide. A perfect correlation would produce a 1:1 ratio (a 45° line). 


Training the network without any input that included S$, =2 attained the output shown 


in Fig. 13. The data with S$, =2 was then input into the trained network. 


3: Transverse Spacing (S,) 


Figure 14 shows that the network could make accurate predictions given 
transverse spacing (S,) different from that with which it was trained. Since there are four 


inputs, the plot displays the output of the neural network versus what the Grimison data 
would provide. A perfect correlation would produce a 1:1 ratio (a 45° line Training the 


network without any input that included S$, =2 attained the output shown in Fig. 14. 


The data with S, = 2 was then input into the trained network. 
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Figure 13: Plot of Neural Network Output (ANN) vs. the Values Attained by the 
Grimison Correlation. 
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Figure 14: Plot of Neural Network Output (ANN) vs. the Values Attained by the 
Grimison Correlation. 
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IV. CORRELATING THE MARTO AND ANDERSON DATA 


A. DESCRIPTION OF EXPERIMENT 


To display the neural network's usefulness, a set of data that had not yet been 
correlated was desired. Boiling data from the Masters of Science Thesis of Anderson 
conducted under the direction of Marto of the Naval Postgraduate School in Monterey, 
CA was not correlated and thus provided an ideal data source. In the study conducted by 
Anderson [Ref. 1], nucleate pool-boiling performance of smooth and finned tubes in pure 
R-114 and R-144/oil mixtures were investigated. The tests were performed with a 
variable number of active tubes (1-5) and variable amounts of oil in the refrigerant (0- 
10% by mass). Note that for input “6,” the number of active tubes does not correspond to 
six active tubes, but rather five active tubes and five dummy pairs. With the apparatus 
(shown in Fig. 15) filled with the R-114/oil mixture, heat was added through the active 
tubes, the simulation heaters and auxiliary heaters. Once saturation conditions were 
reached and the transients had subsided, measurements of temperature and heat flux were 
taken. No correlation of this data was ever developed because of its complexity, although 
interest in its use remained high. 

Figure 15 shows the schematic of the experiment. There are four inputs in this 
experiment: the temperature above saturation A7, the percent oil in the refrigerant, the 
number of active tubes, and whether the tubes are finned or staggered. The output 1s the 


heat flux, g. 
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igure 15: Sectional View of the Evaporator Used in the Marto and Anderson 
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Experiments. Note the Tube Bundle, the Dummy Tubes, and the Heaters [Ref. 1]. 


From first principles, it is understood that as AT increases, g will increase. The 
same is true for whether the tubes are finned or staggered (keeping the rest constant); ¢ 
will increase when using finned tubes. The percent oil in the refrigerant is less exact, 
however. At smaller AZ, increasing the oil % has the effect of increasing g, while at 
larger AJ, it decreases the value of g. Similarly, the increase of active tubes in a bundle 
has a non-linear effect. Vapor-bubble agitation created by lower tubes in a bundle is 
expected to enhance the performance of the upper tubes. The presence of too many 
bubbles may provide insufficient liquid to the surface of the tubes and therefore lower the 


performance [Ref. 7]. 


B. PRESENTATION AND USE OF THE DATA 


Figures 16-18 show plots of the original output data produced by Anderson. In 
Fig. 16, heat flux versus temperature above the saturation level is shown for different 
percent oil by mass. The other parameters (number of active tubes, one tube, and the 
configuration, smooth) are kept constant. Figure 17 shows heat flux versus temperature 
above the saturation level for different number of active tubes. The other parameters 
(percent oil, six percent, and the configuration, smooth) are kept constant. Figure 18 
shows heat flux versus temperature above the saturation level for the different 
configurations. The other parameters (percent oil, zero percent and the number of active 
tubes, one tube) are kept constant. 

If all the combinations of inputs were to be plotted here, another 60 plots would 
be necessary. However, this is not a practical method for using this data. A correlation is 


necessary to make this data useful. Using all the data, the network correlates very well. 


a 


In Tables 4 and 5, the smooth and finned % error is shown. Each value corresponds to 
the average value over the entire group of temperatures over the saturation for each 
configuration. The error reported by Anderson [Ref. 1] 1s given by 


44 _ 0.008 (4.1) 
q 


This is quite an aggressive claim for total error. Even so, however, as shown Tables 4 
and 5, the error of the neural network are well under the error brought about by most 
correlations--not to mention that this ome correlation can be used for 72 different 


configurations with varying temperature. 
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Figure 16: Heat Flux vs. Temperature Above Saturation for Different Percentage 
Oil in the Refrigerant [Ref. 1]. 
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Figure 17: Heat Flux vs. Temperature Above Saturation for Different Numbers of 
Active Tubes in the Bundle [Ref. 1]. 
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Figure 18: Heat Flux vs. 
Configurations [Ref. 1]. 
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Table 4: Average % Error (Over All T) for Every Arrangement of the Smooth 
Configuration. 
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Table 5: Average % Error (Over All T) for Every Arrangement of the Finned 
Configuration. 


9.36 2.69 3.08 2.98 3.48 
2.10 1.37 3.64 | 2.68 


3.21 1.88 
3.96 5.21 2.47 7.20 3.51 





Using the Marto/Anderson data as the target, the neural network 1s able to provide 
accurate results. The network ts trained with one set of data. In order to test its ability to 
predict a system’s behavior, a different set of input 1s given to the network and the output 


of the network 1s compared to other data withheld from the training set. 


l. Temperature Above Saturation 

Data was randomly removed from the set and the network was trained with the 
partial data set. With the input portion of the unused data (the temperature above 
saturation A7, the percent oil in the refrigerant, the number of active tubes, and whether 


the tubes are finned or staggered) was then used as input to the trained network. The 


output portion (the heat flux, g) was used to compare to the output produced by the 
network. The results are shown in Figs. 19 and 20. 
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Figure 19: Plot of the Output of the ANN vs. the Experimentation for Different 
Temperatures Above Saturation. 
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Figure 20: Plot of the Output of the ANN vs. the Experimentation Shown in Detail. 


pas Percent Oil Content 


Figure 21 shows that the network could make accurate predictions given percent 
oil values different from those with which it was trained. Since there are four inputs, the 


plot displays the output of the neural network versus what the Marto/Anderson data 
would provide. A perfect correlation would produce a 1:1 ratio (a 45° line). The output 


shown in Fig. 21 was attained by training the network without any input that included Oil 


%=3. The data with Oil %=3 was then input into the trained network. 
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Figure 21: Plot of the Output of the ANN vs. the Experimentation for Different Oil 
“%o Values. 


3. Number of Active Tubes 


Figure 22 shows that the network could make accurate predictions given number 
of active tube values different from those with which it was trained. Since there are four 


inputs, the plot displays the output of the neural network versus what the Marto/Anderson 
data would provide. A perfect correlation would produce a 1:1 ratio (a 45° line). The 


output shown in Fig. 22 was attained by training the network without any input that 
included Active Tubes=3. The data with Active Tubes=3 was then input into the trained 


network. 
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Figure 22: Plot of the Output of the ANN vs. the Experimentation for Different 
Numbers of Active Tubes. 
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V. CONCLUSIONS & RECOMMENDATIONS 


A. CONCLUSIONS 


Through the correlation of complex data, it has been shown that the neural 
network technique is a viable method for thermal fluid correlations. From the outset of 
this study, the objective has been to develop a better method of correlating data. With 
errors as low as or lower than the estimated error in the data itself, the neural network 
technique displays its potential for complex systems. No longer should 20% error be 
acceptable for correlations. No assumptions were made in the correlation that limits its 
use. This allows for a correlation that can have widespread use with confidence. Finally, 
the correlation involved no "hit-or-miss technique" of development. The weights and 
biases of the neural network are adjusted iteratively giving a correlation without any 


equation. 


B. RECOMMENDATIONS 


1p A more powerful computer could develop extremely accurate correlations 
for even more complex data sets. Larger data sets could be inputted into 
network and a more exact correlation could be found. 

Dy A study could be done on the amount of data necessary to produce 
worthwhile results. Data could be introduced datum by datum. Results of 


the accuracy for each step would provide this information. 
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od 


The area of prediction of thermal fluid parameters can eventually be 
completely understood by neural networks. The area of control can be 
developed for many thermal fluid machines. 

Effort should be made to make the neural network technique more widely 
known and studied at an earlier part of the student's career, as the future of 


heat transfer may very well involve neural networks. 


42 


oP) 


ING, 


ae 


LIST OF REFERENCES 


Anderson, C.L., “Nucleate Pool Boiling Performance of Smooth and Finned Tube 
Bundles in R-113 and R-14/oil Mixtures,” Master's Thesis, Naval Postgraduate 
School, Monterey, California, June 1989. 


Incropera, F.P. and DeWitt, D.P., /undamentals of Heat and Mass Transfer, Fourth 
Ed., John Wiley & Sons, New York, 1996. 


Sen, M. and Yang, K.T., “Applications of Artificial Neural Networks and Genetic 
Algorithms in Thermal Engineering,” Section 4.24, pp. 620-661, The CRC Handbook 
of Thermal Engineering, (editor) Kreith F, 2000. 


Haykin, S., Newral Networks: A Comprehensive Foundation, Second Ed., Prentice 
Hall, Upper Saddle River, New Jersey, 1999. 


Pierson, O. L., “Experimental Investigation of the Influence of Tube Arrangement on 
Convection Heat Transfer and Flow Resistance in Cross Flow of Gases Over Tube 
Banks,” Trans. A.S.M.E. vol. 59, October, 1937, paper PRO-59-6. 


Huge, E.C., “Experimental Investigation of Effects of Equipment Size on Heat 
Transfer and Flow Resistance in Cross-Flow Tube Banks,” Trans. A.S.M.E., vol. 59, 
October, 1937, paper PRO-55-7. 


Demuth, H. and Beale M., “Neural Networks Toolbox: For Use With MATLAB®,” 
Version 3, The Mathworks, Inc., Natick, Ma, 1998. 


Grimison, E.D. “Correlation and Utilation of New Data on Flow Resistance and Heat 
Transfer for Cross Flow of Gases Over Tube Banks,” Trans. A.S.M.E., vol. 59, 
October, 1937, paper PRO-59-8. 


Jambunathan, K., Hartle, S.L., Ashforth-Frost, S. and Fontama, V.N., “Evaluating 
Convective Heat Transfer Coefficients Using Neural Networks”, Jnternational 
Journal of Heat and Mass Transfer, 1996, Vol. 39, No. 11, pp. 2329-2332. 


Lavric, D., Lavric, V. and Woinaroschy, A., “Designing Fin Heat Exchanger With A 
Neural Network,” Revue Roumaine de Chimie, 1995, Vol. 10, No. 2, pp. 561-565. 


Huang, S.-H. and Nelson, R.M., “Delay Time Determination Using An Artificial 
Neural Network,” ASHRAE Transactions, 1994, Vol. 100, No. 1, pp. 831-840. 


Ding, Y and Wong, K.V., “Control Of A Simulated Dual-Temperature Hydronic 


System Using Neural Network Approach,” ASHRAE Transactions, 1990, Vol. 96, 
ative en 7 277 a2. 


43 


[THIS PAGE INTENTIONALLY LEFT BLANK] 


44 


tO 


Ud 


INITIAL DISTRIBUTION LIST 


Defense Technical Information Center......... 00.0.0. o coco e cece ce cece ccc cee cee. 
8725 John J. Kingman Rd., STE 0944 
Ft. Belvoir, Virginia 22060-6218 


Deel ey PROTO tol OUR et sere ites se Paice b sidve vs swalne et pavseas ceedeeseeueececdens ws 
Naval Postgraduate School 

411 Dyer Rd. 

Monterey, California 93943-5101 


Naval/Mechanical Engineering Curricular Office, Code 34 
Naval Postgraduate School 
Monterey, California 93943-5101 


seeeeeeeeeee eee eo eeszseeee ee eae 


oles slit bine yee MCC MCh ty 1.0.5. nanan esses: .2cbseoohe cess sos cc. ses ss ewes 
Mechanical Engineering Department, Code ME/KK 

Naval Postgraduate School 

Monterey, California 93943-5101 


4 voleer ely Beale, oooc do se Rese Men a ee 
University of Notre Dame 

Department of Aerospace and Mechanical Engineering 

Notre Dame, IN 46556 


POTS S50 tank Molen acln CE nt ceo ORIN a Miers a, sands ols casa bale aimee eee 
Department of Aerospace and Mechanical Engineering 

University of Notre Dame 

Notre Dame, IN 46556 


PLOLeSsOs teal MviaU OMT oe ne eM: irish Es 5. vsange vans 044 summa keen eee 
Mechanical Engineering Department, Code ME/MX 

Naval Postgraduate School 

Monterey, California 93943-5101 


JENIN) Wiooyng rte Cree ya lcs alee nee 


6905 North 2200 West Suite 7-F 
Park City, UT 84098 


45 


2 


2 


— 











72 =ae9Ns4 


6/02 22527-200 * 














ong 
ee eee egw Cee « 
ene « 







- ° 
° . 
' or + 
bd *e eeFf “ ) stee 
. “~ a . - 
“e@ « e ~ * 4% 
ay Wane. aa 
ne yom c" «se eg + S elLeee os m 
Pao ae Peat few frees 0 
' SY Pema, 5. aime» s 
° 


. 
- w 
om cas a é 
° 8 Oe We «, ce *. 
e ‘s ~ 0 Me « @ e ww 
a = ied - ° 

of ~ 

+. ~~ ° 


0 We 
=e Nr wy . Ae « ce me 
fe ms ~ w “0 fee ° x ea 
meee a oF coe oe ry 
on ae 


INN 


ae me 
ey 
% 0 wey es “se ta” 
oP et wns cs 
2ee wes 
So meee a Oy, 
MF Pee one 
Seve te, i 1 
owe 5 Ae wie res aw 8. 
Bure % . = re ony Lae ad 
te Oe e. ° s s eA am, « 
ee hl 
Ly = eee 
*. ° 
. a ae iT = . ° . fe 
> e@ we . . -> » 
. A te om =e» 4 ee 
. 
Bree ae . 


ton 
* . e 

. merry ° 
OM © ee ue ; 

e « ch . eee 

Ha mae =e 

oF ations ae 


~. be Se 


I 


. 
. "B00 wncaren oa 
. . * «00, Preys ° wena 
ad ce @ 6: 4 “awl © ae ae 
ae oe Pease . Led Sal 2 ss Cae oo “re tee a, ‘ 
a e. >= o0¢ ® a Cree cote OH hrm os oo @ eeces tes ~ 
tom Sons! 2 ~ an OS @ ee « on, orks - oss “oe CO es oe meee ° ae — 
a Coa fe i = som Ae? 8 -@. 00° ea , 
~ 80 uc om te ~~» Ma ke - man = oot Gon © 0? se ue & Peet Mona Ee a 
. aee ° Cae 2 6 is eee ome etas a eet rid - a ee of SN 0m eee Saray o~nrre of > 
° - €. - se ~ See ~= 
oN tae os °> ee . 
A hi elder Pe & 8% ce secs cae - Py ee S, u sige = be he st Pwamecs POD erenne “ Dye a endee 
: °. 7 - eee hediathed hohe ad teed bd Cie % Saeg Pees 290.8 : mo mes. 
iefostane : Sy stSh ise Veninja a8 Z z : es a. Sea etal Uae sai sace cen So em eeerorees Soxs 
rs Ea ™s = Bo! Oe 8 tee Hoon 
Meee men HA Pens “cong eset we 8 ofe @ Ome asty <2 =s? « - eee . = 
. - = nm 
» >» We . owe re . 
. ado... » > oe se 
” - ° ee 
e ro ° = om op 
e Ostet = Mane , > ame, > © onte 
“se ote- ° foe whe “PrP © wg . 
Oar _— ” ee 
° © © exw s00” e e bad fre 


ee 6 
= se ee to e 
ce, - _—— we mr nm es ° 
e- ana eka =e eee ow | - 
oe = Daw . om *» » ue Ses on 
“* ee =Car ase a sats Dona had 
eee cannes @, acl Po *o%sen o . ewe e Pewerw @ o> te 
fete on OF e 2 om aad ee ee en e ee re ° @ %e 
e Om ane . = =e fe 
~~ . 208 . “2 Od ne 
° tee 
' FAP Fare - 8 s . 
> Pun uous 


- 8e¢ « ce 
“ eile, ae ° 
ee . te 





POF. was 
AAP 2 we, 
bd 3 *Pe%Fes -. ow ax . 
"= ey 7 ere* See . 
Fe em * 0 Pome Aus ae 
4°2™*, Stee *"ay & 


DUDLEY KNOX LIBRARY 


WN 


4 ng ote * PnOrgue 
AW *® eae ot ~~. os 
we Rad fa MODE Oe Oe oy On Wod.s _ 5 eutee 
Pa~p odr ene 8 La en, 5 heh ta ee te WO te a0 2. 
= © wuss NTRP a8. woes, thet SE ee and $a ae weer, ot, | 
a a = ol ow bs cewanre . 
- we, 


o com =e 
° 2 @ OL PS®r.w « . * 
“Ma © a6 .0%5 oo * ot" De. 
. Fy Ee re Sy Fae ore 
oS #8 are eteu f Levee etme vc wen, °. == nngee 
ee rat » 
Li © 0c, = = = Comm ° fo we 
eco 8 ve ee nt we .? oo «- 
oP b= veage ee om . . © * eas 
= @ee oe & ee - e se bd 
. . . * sad ane . ° ee . 
% =? oP ° & oPecs a Lt o © Brag 
—_ . - e. 8m e. = = 6 
ee ee, e 


oom OF Loe 


3 2768 00403388 6 


° "ae 
eae - men 

sP St oe Oe on 

“e OP Psa, 
. ee . 

OF us owns 





~~ 


90 w-siwe ay zs “SOP BC wad ny a 
“te C.me *. = - on Pe Sealed -detake ee tO aind pets J toast sep ppt ee ae ete 
Petey 6 bale « i comes SO Cas om a —— i & wren Oe ae oe vee qecr~e 
CRP worm °# ° 86m see o e%e oe ome on. MO) oes ose ae ae a a BO SOE @ 
Sot ad . ey on? S250 ea 2, ts 200 . =, . = 2 Se Gorm 
A eS ears 6% OMS eke, = Sa saeteae cs onus a were ok: ve 
— PH letetianw, ete <a. ey © Meta, 2. 
e+ me % shen, Pre Coca "eet xm, 2 Mery wwe, rete 
ed ete ~ se © Poate 4a. ee @ 2 om, 2? « “om th 
Ps . - . ® ee Row * were rR ame ae 
—_ * wee = ° Boe our @ecewy 
« . =. 2 e= was pebn wony @ Tehran. « te ee “f 
' e o> Sesernin os Oe ews we one - SiGe  sivtsese « 8 
”~ ny Poms iw ° to 08 “ewe 8 105, Pen o ~ . provvene a 
FP” Pose « . 00 Ae ecu ene ors ' Oe we "8 »® i Stews 6 Oke ma ww &, . @ pease ‘ es ww me 
ems, » 292 oe we | %es “ o s = ‘8 * e e » a2 & Bote Mise o S*dences 

vee » . “ . » wee Fs aie 
20 a ay ® aie 455 pute - > bee «@ & mean 

¢ ema e oe 8 ta mt 7 . oun? » ose 7 om 
cote bay aCe ao - = 6 e ° " wa . at me ay ° 
Wee Poh e cons = Fern ON @ ae le 2S ke Pe eee Sterme 6k fk . 
. ry Bee @> een ee 
* @e ce «e 


. 
ov. mem on 
Foor 3 oe » * s . Aste & 
ww 0 PsP ete. sea. S ° oe o. wa : - - _. a 
e et te > ane & rem bed ees en . e 
. te Pa Ree = eae 
& Seren ony be 7 cage = Sin 
ph OS ete tare aa ~. 
i a oe . . ° s 
== 3.6 


e 
° 


ry > fe 2 
- & 


eee 
» oe. . 

F8mG? 
"e aOe . ates. + Lhd 
° 

eee ew Pate © et) 
©¢ @ewe 
© tees 


Brew « ° ¢ weerse eran a ae 
le Ty oo? @s%; ae @erte, - tf e ® 
eae We & ver ay Wen" we ®» * e@ «0 mo 
ee ee e . ene og Ve ee s 
ory a Bre 8 erwreton aise = onth 
baited me OPN oP eri aes “s. eon ee ow 
o 0 lw, ad) et > © %e@ =. ~e 
~ ure os o, 


2° ee 
- 9 0 vtrmm 15 = 
P tee pment oe O- Rom Oe teneg on ~ 
iS Oe" ae ow ey) ee a OOS whe we.e un - ewe 
O omer (Pees Aan, Or. wm o Mw oy ®et— 0. 
i Swe we Femeet F*oteagt = ° 
eal conic, Lee ce SN ow O08 - oS SH o-temmeee 
OER OP OF aoe whee Ne car oe 
. he-dae L e 3” «© oem urate maa « es 

. tee ° %% me ty me POPE cageer a 
= 4s Scere, & % = - ve eo ee e . 
7° e ehese A ae ate 


2 os e. 
© Se wert oe < ges 
Weer eo « POe cy 
. We 6 ° Peeee » s we 
ee . ee ue = a "9% be 6, 
» od 2 AS teres wie Teh pree bd he Os 
%. = © es eee eft ce © Set mere, ar i . 6 rE oe pw wucata a 
e er + @.. oy ° @ o ww ° ec 8 wee i ae | eee. etm. o we Sow eas « 
° s 2 c°. 8 ee RNP erm oP ° * «© Se tePS ers, Mes “gece "OPM ot tne ny ois t we, te mane cers 8° 0 th © meee : 
See) Oe) Seana ee es % wee oe oP ree Cr 4 ote fe Page, ee > ae, 0 MAT e antag ‘ 
ot wo Fe ewe me os welt Ce of SF 6 toh og woe shy ‘. en bet Oe wee, 7 ER A . ° 
ee 4 *e ° . » © 0 ® a.eu, eos 8 ®% @ .. . Xe es Sie ate ee mao s et etete se, 
Scare sre & oe Ss. ae oe esis Ret ee tl oes coy we ON Ne a. om ~ ce Ma ee =, BIS Owe chee. < ofa eg ; % , 
Pee bares wer =e te poe POP we of iw 9p © $25 hts - aa . we elew + * eters aye o = «~ ote of . 
See tee, -, tw om my @ Pus Fe ty o wae feo wee t é owe ee eo oe, are e Cr men ’ 
Pete RA ee bebateaeh (i>) uae ° © oot Fewer, 2s Ee og ° & ee. eo ae *@e =s 98, oy w es 20m neeem gus weSeub ewe 2 
el SSM, ee > = 660 © sts * 0) eepme. oy etme a.e,e ” ° 4 © des 
"Wee ee ween Cote ron, rh: ror . 5 Mew Sees oaay es we at ) > . wok -~ ee 3 + Memeo e eens es = 
oP mee wy . Fees oe e © @rre wae Os 8 Ow 2 Steg i ° 2 oe, z ee 9% —. 
© @ Weer o © Me ~ en. we 6 ie ee vs ae r nm ° @ . 
-_ «© mee er wee we og Som ot”, os . oe eee o othe e 
70 comes. fy ° ° or See oe * 
¢ On ee 
sen Fh ymek es Cw ty © es ow Oo oe we 
= a te wee Mp ae SS 
pe ee ye 


= © ate. 
. F e Per» O.e Oe yt ways "amet 
7 wees 0 MIs « 3 
* 28 O08 ome Oe Seem a 
te te? @ Oy _ “ peal 
ewes 20M ome y) e we oy ome wo s i rs rh oe OO + teen eps Te. ry etek ae, 
x mm ees, “oS tee a cme =. . netheg” 5 Om wewres see Woe 9-2 an, 

peta 0 Ce os ° we © neo) Nts we eS “ort wegen, Set Neuse 
’ - sexi e fit © se bh a © ae? cwwe 8 a) . 
: my oe “wee ye ae ° fre ae evae z ose ~ ee 2 4s ss teens 
ary "Ye ~ 06 os o> © rh oe PU ae 8, om -¢ © ow 
= a oe ° « &e Lets OOo om met ms g pone om Bierce — Pe. 
—=e-. A) ° =e e8 PO Faas e e %e 

% se ®e "s20@ ®% te 84 He © 9 P wre 
One!) 2O ra es ‘o Ow bd od - e 28 Te fee 
° e ee ° & te be .gapeste © 8 eee tes ~*~, é 
. coe oo etry. ea weg © e bad} - 

w Pe ° 

es . 


° ny 
i ° 2 <e ong oP ae Bee 
Lo: a F “ 
wae? ee 2 ges qos Tt 
ee 
se . Pea, 
see 


e0 we °* Nes 
O° 8e teen, orca bai * e8 
omen. 


tad 
-- 
© & etems ort 





ofa, m0, 
ase ” 


ete 
"et 


ee em s 
e™ @ 6 ote elds id . 


Te m& aoe 
Ce a ee « 
— @ ce 


2@. Sedge ame dy: 
bind yh bn alt MO 
oheneges Po ereh 





ee « our rw bad e 
Or & oe vas seo 0 bet?) we oe Wars 
°. = we eat ee eT oe. oo 
Rtoroe® vp Oe eters Pane wate ie ° & e,, 
ee te’ Seb ww « 
Me Cons vo, wag 

Wes oO © oe 

: om OPE wae cures se 
? oO 4 ata rrmerys Dcocep : 
Ss = ctor beam, pddaseabe dt) ra soe yes ey 
FOF Oho we me fem Co) res teow ad uy 
owe © Ate OUP enh tet cco e 


=. ae “et 
Belial PLOY Hn ererarrant ST .eeP bes 


e se 
crete eatin SSeS 
A at od 7 Em oe way SMetepech oe ~ Lt rapa 
-! o- eae “a Lr 7 Py 
i C FO em ss * - Fran” oo) Btbllsewon neces Nee 4 ep © ates Uwe 
ee 68 Seem os ee ce tA Bee. CL Pra fa- mt atge or 2 
° i ° @ ° 6 we see te ee ’ Hae > oe@ tore © te tg « 

08 c08 e Ctre 4 yl e 6 * ee ce ote cee Seth BS Beye = a & . . “a 

° woe etd oo, Mh ee one aa ae SOne moar, | ao oe Re pe Pee oe io re = wee 
a Fe Peer, Pre ro ake ©o wet hesty ACA ne > e vw . ode Neo se ge Tite tone 
ee dB oe wt. bd Mle ¢ . %e a eee eo 3 a a 

© be on. ew Arena a, = Wy we oy ee i . “ure tee, 
= oP & oa ee SO whee we ° -teo ow =f > se mee > + = 
+ "Fe mame Se sey “tots om, ‘ot erree &e fete We *.- @ oe te 
Ory "VWeliw, fep, =e cme =e we . ee we. . © Oo oe ete tees + . m= . 
Peryorny, wolereee, oben ys eo wen. ge , - us = nw . b . ° ee * Siete 
ioe | Yu. as ohms @ © © 09 B® oy Pe, . ee e ; - “= we woe w wbia.e e 
% Gee Oem 0 cre), e ees. Ba - te 
2 © 88s ost «. ve, 
* Ose, s0r. ,, 
"Peuterk’ ia, 


WUE Bee, 
: Srownca ce cote. 
ee Be ® a5 %, 

ove 

8° Peer sig se * SP: a 

© 6 ons Fouts, @ ee 

to ot%e * 

* Lad ce 

- ew% ae te Fe 

*@ meer © ween 

oo at) ae 

®-as ». | ny 14 

@ mn 8 tee wr 

Werte oe %.° e 

were 2. e ©Be 004° Qeer 

“Wie Bf 0a. Gc0ea. Bee 


oh. 
stem mage ° Sea =e, 
e *E=~ «¢ 
tb Ao e : ks 
Sere he 
© Mwe seeus at 
e = SMe een 
"28 @ . we . 
ee eo 6 3 . 
Resco se O° FP o%e.wetes « ou 
Bee & eee 
“Se. 


o8. feet oer cnt 
2 Ws Sorat tye 
"Cee Ce 
we 





- 
oo Wee vise 


o_ 
oer koe 
° - = <a le oe 
7 Ot Pom Ui re To eee 
ad ° ® ° 
e ie e.0@ .e .. 

theo . by igage 

Eee “ty cee foe, 00 oP Pans. 







= bow, J 

ate ibe OL aed Se 1 6e2e ope ore ore. teq o,, 

2M FS Rte wesee ° oo wBene “ *Moeme s te 

oe We Pe 2009 05: «war nn * Se wyerey oo tthoctatea’ ep Appel 

ee © © CA -Per.esipcarre ®-e° ere. Fw Rhee ta TTT 

wee . Vee waren Soe, wR Od oo 0 
, 9 Gone 


death tl Ty 
Meee, 


Ss MOP Mehes. 
Bon pePone, o°seg 





e ° bd ® 
= e 
2 t.- = 
> © .e © 0% oben seep = a ews ¢ » 
2 ee 8, A MGs 5 os bam "e 
Le) ° @ os ™~ a *4t *. Mee e 
Seo 5 Ore, te ® 
0 @@e ay ms 4s bs 

. ied Oe" were ores 

. O44" Ge . b +e 08 
= Os re tah og Pare 
fc int To) 8 OP %0 “erste be ange " 


e 
ee J 
ee e id « 
A te *° =e. = 
Qu? ee ee ti ® s »-> 6 
- ee - *e te ee ee ower 
wy et) Le sk at) se = 4 e e ‘ 
*e ae “ ® « 
ver @ seg ; . o as o- we 8 2 - ® ry 
es ee ac a) ° ° ee ee ° ee 
_ 40% © © = & et aio "Oe sae 5 . 2 Peor 2° os ° 
° .  .wteec tg — er r ° 
t CeD amas ‘. a Petes ef » we 0e* te, » @e eo « we ras ° i 
no mage Hon a 7 “ se . O° wore, eh . Ole eae if e e e 
Rat btn de 1d Oo eee ow Cee owe s- =o ek * nee #0 es > ee Oran e. 0 et « ee . eas 
pea ree coe yee cs, eee oot arn Bare, y 5 . Secs wee esac Powe pes aw vera cae oe Se 0:81 Ae 
Soe ae ener ese ap edeare »: 60° Cet! ale Gem. , Ess Oy . Samy, "et ee eee 
se wer ne? Rite ft Gl eats , <e 
oe = Oe ON cw PS thanehed be BCT ep ime 
we mee erent 





© Pee ee. 
5) *Biseens an att 
”~ ad, ee 
PW Pl cet ngtges 1), OM? caeha Rh, cota 


es ee @ © a. 
2 2 

> _ eo 

sd wees = 
a 


